The angular momentum of a magnetically trapped atomic condensate 
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For an atomic condensate in an axially symmetric magnetic trap, the sum of the axial components 
of the orbital angular momentum and the hyperfine spin is conserved. Inside an loffe-Pritchard trap 
(IPT) whose magnetic field (B-field) is not axially symmetric, the difference of the two becomes 
surprisingly conserved. In this paper we investigate the relationship between the values of the sum 
or difference angular momentums for an atomic condensate inside a magnetic trap and the associated 
gauge potential induced by the adiabatic approximation. Our result provides significant new insight 
into the vorticity of magnetically trapped atomic quantum gases. 
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A magnetic trap constitutes one of the key enabling 
technologies for the recent successes in atomic quan- 
tum gases 1]. The most commonly employed mag- 
netic traps includes the quadrupole trap (QT) as in the 
magneto-optical-trap (MOT) configuration 2] and the 
loffe-Pritchard trap (IPT) 3]. The direction of the B- 
field forming the magnetic trap is generally a function of 
the spatial position. For a trapped atom, its hyperfine 
spin adiabatically follows the changing B-ficld direction, 
and the atom remains aligned (or anti-aligned) with re- 
spect to the local B-field. As a result of the adiabatic 
approximation, the center of mass motion of a magneti- 
cally trapped atom experiences an induced gauge poten- 
tial from the changing B-field [3, 0] . 

In a variety of magnetic traps, e.g., in a QT, the B-field 
is invariant with respect to rotations along a fixed z-axis. 
As the cse of single particle dynamics [6], such a SO (2) 
symmetry leads to the conservation of the ^-component 
Jz (= Lz + F^) of the total angular momentum or the 
sum of the z-components of the atomic spatial angular 
momentum L and the hyperfine spin F. A different sym- 
metry exists for an IPT giving rise to a corresponding 
conserved quantity (= — F^), the difference of 
Lz and F^. To our knowledge, this surprising property 
has never before been identified explicitly. We thus feel 
obliged to present this letter because it significantly af- 
fects the vortical properties of a global condensate ground 
state in a magnetic trap. 

Our work is focused on a detailed investigation of the 
relationship between the gauge potential and the associ- 
ated values of Jz {Dz) in a magnetic trap. For a spin-_F 
condensate, due to the appearance of the adiabatic gauge 
potential, the possible values of Jz or Dz are restricted to 
a definite region [—F, F]. The gauge potential of our for- 
mulation is directly related to the effective trap rotation 
studied earlier in Ref. While Ho and Shenoy mainly 
studied the orbital angular momentum component in an 
IPT instead we concentrate on the conserved quan- 
tity, the sum (Jz) or difference (Dz) for a QT or a IPT. 



This paper is organized as follows. We first consider 
a spin-1 condensate in a QT. Making use of an effective 
energy functional appropriate for the adiabatic approxi- 
mation 01, we prove Jz G [—1, 1] with the actual value 
determined by the angle between the z-axis and the di- 
rection of the B-field. We then generalize to the spin-i^ 
case. Finally, our result is extended to an IPT. 

The Hamiltonian of a spin-1 atomic condensate with 
N atoms in a magnetic trap is H = Hs -\- Hi with the 
single atom part 



V>^(f) 



-1^ + l^BgpBirjF ■h{r) 



ip{r)dr, 



and the atom-atom interaction Hamiltonian 



Hr 



ipir) — (r), '00^'' (r^, V'|;^''(r^]'^ denotes the annihila- 
tion field operator for the z-quantized F^-component of 
m,n,p,q = 0,±1. M is the atomic mass, and /ib is the 
Bohr magneton. B{r) and n{r) denote the strength and 
direction of the local B-field. ft = 1 is assumed. The 
Lande g factor is 5f=i = ~l/2- 

Within the mean field approximation, the field opera- 
tor 'il'{f) is replaced by its average {'ip{r)). To introduce 
the adiabatic approximation, we define a group of nor- 
malized scalar wave functions (^„(r): 



(1) 



6=0, ±1 



where £,^{b,'r) is the eigenstate of the i3-quantized spin 
component F ■ n{f) with eigenvalue b, satisfying the rela- 
tions F-n(f) ^^(6, f) = b^^{b,r^ and ^^Hb,r)^^ {b' = 
db b' ■ In the z-quantizcd representation, it takes the form 
[b, rO = K-i (b, f) , (6, f) , 4^ (6, f)]T . In this study, 
it is important to distinguish (&, r) from the eigen- 
states ^^(0, ±1) of Fz with eigenvalues 0, ±1. In explicit 



2 



form, we have ^^(-1) = [1,0,0]'^, Cz(0) = [0,1,0]'^, and 
C.(1) = [0,0,1]T. 

A magnetic dipole precesses around the direction of a 
B-ficld. Majorona transitions between different ^^{bjr) 
states can be neglected when B{r) is large enough. Thus 
the atomic hyperfine spin adiabatically freezes in the low- 
field seeking state £^^{—l,f) during the trapped center 
of mass motion, and ip-i(r) — if(r) and ipQ^+i{r) — 0. 
Similarly the z-quantized mean field becomes 

i:{r) « (V3(f)) - V]VC^(-l,r)^(r), (2) 

with lyj(r) a B-quantized scalar function. 

Substituting Eq. ^ into the expression of H, we can 
obtain the expression of the condensate energy Eg_^ as 
a functional of the scalar wave function '■fiifY [ip\ = 
N J if* {r)£ad^{r)df. Here fad defined as [J, 

gad ^ ^-'""^.f^^' + + + K + f IV^P. (3) 

The gauge potential induced by the adiabatic approx- 
imation is A{r) = -i^^t(_i^^v^^(-l,7=) and W = 
{\V^^''i-l,f)V^^i-l,f)\-A-A)/{2M). The trap poten- 
tial HBB{r) is augmented by Vo{p, z) from other sources, 
e.g., an optical potential with rotational symmetry along 
the z-axis. Under the approximation of contact pseudo- 
potentials, trapped atoms collide in the same spin aligned 
state. Thus the 172 term is proportional to the scattering 
length 02 in the total spin channel i*tot = 2. 

We show below that the gauge potential A{r) con- 
strains the values of Jz for a condensate ground state. 
In cylindrical coordinate (p, (/), z), the B-field of a QT is 
expressed as B{r) = B'{pep — 2zez)- An optical poten- 
tial Vo{p,z) is introduced to push atoms away from the 
region of small B{f), eliminating the deadly Majorona 
transitions 0, Q. Because this B-field is cylindrically 
symmetric, F ■ h{r) commutes with the z-component of 
the total atomic angular momentum Jz — Lz + Fz- This 
allows us to choose £^^{~l,r) as the common eigenstate 
of F ■ h{r) and Jz with any possible eigenvalues. In this 
study we choose ^^(— l,r) for convenience to satisfy 

Jze{-l,r^ = ^- (4) 

This constraint on Jz limits the respective values for Lz 
or Fz, in a sense equivalent to a gauge choice. The ac- 
tual value of Jz for the ground state of a condensate is 
determined by appropriate system parameters. Explic- 
itly, a simple rotation gives ^^(—1, r) = exp[z0] exp[— iF- 
e^P{p, z)]^z{—^)- The angle between the z-axis and n{f) 
is P{p,z) = arccos[— 2z/ \/ p"^ + iz^]. 

The _B-quantized ground state scalar mean field wave 
function is denoted as fg{r), determined from a mini- 
mization of i?ad[</5]- The cylindrical symmetry of both 
B{r) and Vo assures E^^difip, <l>, z)] = E^dlifiip, 4' + S, z)] 
for any 9. Due to this SO (2) symmetry, if there exists 



only one normalized ground state (up to global phase fac- 
tors) as in the situation considered here for a condensate, 
it has to be a common eigenstate for all rotation opera- 
tors exp[9d/{d(j))], i.e., an eigenstate of id/{d(j>). Thus, 
we take Lpg{'r) = (pg{p, z) eiq){is(l)). On substituting into 
Eq. (21), we obtain the z-quantized mean field {'4>{r))g for 
the ground state 

{i^(r))g - ViVe^(-l,rO^<,(p,z)exp(z5(/.), (5) 

which is an eigenstate of Jz with an eigenvalue s, i.e., 

jS{r))g = s{Mr))g, (6) 

because of Eq. (jH). 

Equation ([5]) and the expansion {i^{r))g = 

Em=o,±i(^™^(^)s^^('") gives 

(V'i?>9 = V7Ve^(-l,r^^,(r) = 6„(p,z)^,(p,z)e^(^-™)^(7) 

with b,n = ^l{m)e-'^y'^''P-^^£,z{~l). The individual spin 
components of state Eq. ([7]) naturally carry topological 
windings as a direct result of the conservation of Jz . 

To determine the value of s or Jz for the ground 
state, wc first compute the gauge potential A{r) for a 
QT. With the expression of ^^(— l,r), we find A{r') = 
cos I3{p, z)e^/ p. Using the expressions for A{r), fad, and 
i?ad, it is easy to show that for a scalar wave function 
(pg{p, z) expiiirKp) with any integer m, the energy func- 
tional i?ad satisfies 

Fad [^g (P, 2)e*™"^] = Fad [^g] + AF„ [f>g] , (8) 

withF,„[<^g] = / df1(^;(p,z)|2[(m2-H4mcos/3)/(2Mp2)]. 
In addition to the centrifugal term proportional to , a 
term linear in m appears due to the A ■ V term in fad- 

In the following we show that the above linear term 
is important for the value s, which we determine with 
a variational approach. Because Fad takes its minimal 
value in the state (fgexp[is(j)], we have Es.di'fge^^'^] < 
-Ead['^5e*("±i)'^]. Together with Eq. we find the nec- 
essary condition satisfied by s: AEs[ipg] < AEs±i['f>g] or 
|s -I- C| < 1/2, where the coefficient C is defined as 

/dr1^;(p,z)Hcos/3(p,z)/p^] 

Jdr\0;{p,zW[i/p^] ■ 

When the correlation between cos/3(p, z) and p^^ 
is neglected, the factor is approximated by C ~ 
J (ir|(^*(/9, z)p cos/3(p, z), i.e., the expectation value of 
cos f3{p,z) for state (pg. Since |C| < 1, we have |s + C| > 
\s\ — 1, which is the same as |s| < 1 or s e [—1,1]. 
Without the gauge potential A{r), we would always have 
AF„J(^g] = J dr\if*g{p,z)\^{2Mp^)-^ > 0. Thus the 
value of s definitely would be zero. Therefore, the ap- 
pearance of a non-zero valued s G [^1,1] arises due to 
the induced gauge potential. 
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We now generalize our result to atoms with an arbi- 
trary F and inside any axially symmetric B-fields. Anal- 
ogously we can prove that the value s of Jj in the ground 
state satisfies the necessary condition 

\s-7^fFC\<1/2, (9) 

and s E [~F,F] with rjp — sign{gF). The result of s e 
[—F, F] and the conservation of is independent of the 
form of the atomic interaction potential. Although its 
strength g2 docs affect the wave function shape, thus can 
influence the value of s through the factor C. 

The condition Eq. ([9]) also allows for a rough esti- 
mate of L^. A straightforward calculation gives (L^) — 
s — rjpF J dr\(p*{p, z)\'^ cos l3{p, z) for the spinor mean 

field {tp{r))g. Neglecting the correlation between 
and cos I3{p,z) as before, the value of (L^) becomes 
approximately s — tjfFC, which lies always in the re- 
gion [—1/2,1/2] according to Eq. ([9]). Therefore, the 
value (Lz), or the weighted average of the winding num- 
bers, is generally a small number, despite the wind- 
ing number s — m itself, for the component {ipm {r))g, 
may take any integer in the region [—2F,2F]. We find 
(Fz) = rjpF J dr\ip*g{p, z)]"^ cos P{p, z) from the expres- 
sion of (Lz), a qualitative reflection that atomic hyper- 
fine spin is ahgned {gp > 0) or anti-ahgned {gp < 0) 
with respect to the local B-field. 

Our result above allows for the direct creation of vortex 
states in a quadrupole trapped atomic condensate. For 
example, assume a spin-1 condensate in a QT plus an 
"optical plug" Q satisfies Vo{p,z) — Vo{p,—z), then we 
find C — and s — due to the spatial reflection symme- 
try about the x — y plane. The ground state components 
(V'±i(^))g then automatically carry persistent currents 
with winding numbers ^1 according to Eq. ([7]). In addi- 
tion, the low field seeking atoms are trapped near the x-y 
plane at z = because 1 5(^0 1 is an increasing function 
of z. The populations for the three z-quantized states, 
determined by r) and 1, r), are of the same 

order of magnitudes. Therefore, when a ground state 
condensate in the "plugged" QT is created, its ±1 com- 

" ( z^ 

ponents {ipj.-[{f))g are single quantized vortex states and 
can be directly resolved with a Stern-Gerlach B-field as 
used in Ref. 

The qualitative example above is confirmed by the nu- 
merical solution for a condensate of 5 x 10^ ^^Na atoms 
in a quadrupole plus a plug trap. We take B' = 22 
Gauss/cm and Vo — J7o exp[— p^/cr^] with Uo — (27r)8 x 
10'*Hz and a = TApm. The ground state distribution 
p, = / dr\{'ip'l^\r))\^ is found to be p±i = 27.2% and 
Po = 45.6%. The phase and density distributions for the 

I z ) 

three components {ipQ are shown in Fig. [l](a, b). 

We also can expand the ground state {'4'ir))g in terms 
of the eigenstates ^xi'/n) of F^ with eigenvalues to: 
= E,n=o,±i^i^^rn\r))gUm). We then im- 
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FIG. 1: (Color online), (a) The phases of the z-quantized 

components (i/ji^') (left), (V'o^'') (middle), and (right 
panel), as functions of the azimuth angle 0; (b) the density 
distributions \{tij["^)f (left), K^/i^"')^ (middle), and 
(right panel) of the z-quantized components as functions of 
p and z; (c) the integrated density distributions f \{il)[^^)\'^dz 
(left), J\{i,l''^)\'^dz (middle), and f\{i>'Ll)fdz (right panel) 
of the a;-quantized components as functions of x and y. The 
units for x, y, p, and z in (b) and (c) are all arbitrary. 

mediately note that {il>m\r))g is a superposition of vor- 
tex states with definite winding numbers or ±1, e.g., 
0o"\r))g = {VN/V2)[b-i{p,z)e'-i' -bi{p,z)e-% The 
density distribution of (■0o'^±i(^))g as shown in Fig. [Tfc) 
clearly illustrates the interference pattern along the 
direction. As is demonstrated in Fig. [T]^c), the mid- 
dle panel for \{tpa'\'r))\'^ clearly displays the double peak 
structure along the azimuthal direction, arising from the 
interference of the terms proportional to e*^"^. Thus, if a 
Stern-Gerlach B-field is used to separate the x-quantized 
components {ipm\r})g, a superposition of vortices with 
different winding numbers would be obtained. 

We now extend our result for an axially symmetric 
magnetic trap to the widely used IPT whose B-field pos- 
sesses a different symmetry. In the region near the z- 
axis, B{'r} = B'[cos{2(j))ep — sm{2(j>)e^ + hcz], the an- 
gle (3(p, z) between the local B-field and the z-axis sat- 
isfies cos/3(p, z) = hj \J p^ ^h? . In this case is no 
longer conserved due to the lack of the SO (2) symme- 
try. However, we find that Dz is now conserved be- 
cause it commutes with F ■ B(r). Therefore, we can se- 
lect the low field seeking hyperfine spin state ^^{rjpF, r) 
as the eigenstate of Dz with an eigenvalue —r]pF, the 
same spin state as used in [Hj, again defined through 
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a rotation {i]FF,r) = exp[—iF ■ fi±P]£,z{r]FF). For 
/i > 0, we find the induced gauge potential becomes 
A{f) = -ripFil - cos P{p, z))e^/p. 

Adopting the same notation as before we denote 
{'i>ij))g — "/Nipg{r)S,^ irjFF.r). Interestingly we find 
Ea.d['~p{Pi 4>, z)] — Ei^ix[ip{p, (p+O, z)] remains satisfied, and 
the ground state takes the form ipg = (pg{p, z) exp{iu(p). 
Therefore {i^(r))g is the eigenstate of with an eigen- 
value d = u — rjpF, and its components {'ipm\r))g = 
•\/]V6^(/ci, z)e*'™"'""~'"'^''^ analogously carry a persistent 
current with a winding number m + u ~ "qpF. Here 
h'^ = ^|(m)e~*^f''^z(77i?F). This result is consistent with 
the ground state vortex phase diagram for an = 1 con- 
densate found numerically in the z = plane of an IPT. 
The conservation of as found by us, however, calls for 
a simpler labelling of each vortex phase because only one 
of three integers (m i , mo, 7n_i) is independent, as with 
Eq. (15) of Ref. [13]. 

Following the same reasoning as before, we find 



\d + riFF{l - C)\ < 1/2, 



(10) 



for d ^ tjfF, and d g [—F, F] or the value of in the 
ground state Hes in the region [— F, F]. 

In an IPT, atoms are trapped near the z-axis where 
the B-field is essentially along the z-axis direction and 
^^{r]FF,r) is approximately the eigenstate ^^{tjfF). Lz 
then is essentially always zero corresponding to a ground 
state without a vortex. The angular momentum differ- 
ence Dz then becomes d = —tjfF. 



Several previous proposals 11| and experiments |12l | 
on creating vortex states unknowingly have used the fact 
that (V'(r))g in an IPT is an eigenstate of Dz- For exam- 
ple, in the experiment of Ref. [12l |. spin-1 atoms initially 
were prepared in the internal state ^^(—1) {Fz = —1) 
with no vorticity in its spatial mode. This corresponds 
to Dz = 1. To create a vortex state, the bias field along 
the z-axis was adiabatically inverted from the +z to the 
— z direction. In this process the atomic internal state 
was changed from ^^(—1) to ^^(1). If the whole opera- 
tion is adiabatic, the commutator [F-i?(r, t),Dz] = will 
be maintained. Consequently Dz is conserved. In the end 
when the internal state was changed to C^(l) (Fz = 1), 
its orbital angular momentum became Lz = Dz + F\ — 2 
or a double vortex spatial mode as was observed |12| . 

In Ref. when the bias field is adiabatically switched 
off, the direction of the B-field adiabatically changes to 
lie in the x-y plane. Dz is again conserved during this 
process (= 1). When complete, the atomic hyperfine 
spin state is changed from C(— 1) to a superposition of 
all three components ^^(0,±1). As a result, different Fz 
components are then associated with vortex states with 
corresponding Lz — Dz -f Fz- When the three internal 
states are separated, two of them are observed to contain 
vortices with non-zero winding numbers. 

Creating a ground state condensate with Dz 7^ —ijfF 



is quite challenging inside an IPT. This was considered 
quantitatively by Ho and Shenoy who obtained an 
approximate gauge potential resembling an effective trap 
rotation when expanded to the first order of p. How- 
ever, their expansion easily fails away from the z-axis. 
A numerical discussion on this challenge is provided in 
(lo| . The more general constraint of Dz ^ —i^fF, i. e., 
the necessary condition Eq. (llOp of the angular momen- 
tum difference Dz, was not obtained in p, From 
a direct calculation, it can be proved that if C is ap- 
proximated by J dr\ip*{p, z)\'^ cos (3 (p, z) and d 7^ ripF, 
we have (Lz) — d + ijfF J df|i^*(p, z)p[l — cos/3(p, z)], 
which can never be greater than 1/2 according to Eq. 

m. 



In summary, we have investigated the angular momen- 
tum of a magnetically trapped condensate. Inside an 
axially symmetric trap such as a QT, the total angular 
momentum Jz along the symmetric z-axis is found to be 
conserved, while the angular momentum difference Dz is 
conserved in an IPT. Both conservation laws reflect the 
underlying symmetries of the traps' magnetic flelds, and 
the values of Jz or Dz in the ground states are determined 
by the gauge potential A{r). In the global ground state, 
the corresponding eigenvalues of Jz and Dz are limited 
to G [—F, F] with the precise values directly related to 
the angle between the local B-field and the z-axis. 

Our results provide significant insights into the study 
of magnetically trapped condensates. The conservation 
laws we discuss reveal an important observational con- 
sequence: the components {ipm {'r))g of a condensate 
ground state {ipi'>^)g automatically carry persistent cur- 
rents with different winding numbers. Furthermore, ac- 
cording to the conditions Eqs. ([9]) and (fTO)l . the values 
of Jz or Dz, or the topological winding numbers of the 
spatial wave functions, can be controlled through the an- 
gle /3{p, z) between the local B-field and the z-axis. We 
have shown explicitly for a condensate in a QT with an 
optical plug, where the atomic populations for the 2F+1 
components have approximately the same order of mag- 
nitude. Therefore, vortex states can be present already 
in a ground state condensate without requiring adiabatic 
operations as in [3, [l2| . 
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